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Abstract 

We prove the universality of correlation functions of chiral unitary and unitary ensembles of ran- 
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local asymptotics around the origin. The possible physical interpretation as the universality of 
soft spectrum of the Dirac operator is briefly discussed. 
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1 Introduction 



Consider a four-dimensional gauge theory coupled to fermions in given representations of the 
gauge group, and with given global (and possibly spontaneously broken) flavor symmetries. In 
the chiral limit, the associated Dirac operator \]p anticommutes with 75, 

{i#7s} = . (1.1) 

This means that the associated eigenvalues of ifl, defined by ilp^> n = ^n<fin appear symmetrically 
around zero: ±A n . The corresponding eigenfunctions are, for X n ^ 0, (f) n and 75 (f) n . The accu- 
mulation of eigenvalues near A = determines whether or not chiral symmetry is spontaneously 
broken through the formation of a chiral condensate O. The study of the Dirac operator spec- 
trum near A = is thus of fundamental importance for our understanding of chiral symmetry 
breaking in gauge theories. For example, for QCD with Nf flavors it is expected to lead to the 
conventional SV(Nf)-^ x SU(./V^)r, — ► SV(Nf)y scenario, with all its implications in terms of 
low-energy effective Lagrangians. 

In a series of papers [||, ||, f|, Verbaarschot and collaborators have added substantial new insight 
to this issue. Their central assertion is that the spectral density of the Dirac operator very close 
to the origin A = should be universal, depending only on the symmetries in question. One 
startling consequence of this conjecture is that the spectral density of the Dirac operator near the 
origin need not be computed in the gauge theory at all, but can be extracted from much simpler 
random matrix theories. The pertinent random matrix ensemble is determined by symmetry 
arguments alone. To find universal features of the Dirac spectrum near the origin, it is essential 
to first consider the problem in a finite volume V4, corresponding in matrix- model terms to finite- 
size matrices, and then magnify the spectrum near the origin at a scale of order I/V4. The usual 
spectral density p(X) is defined by the following average over gauge field configurations: 

p(A) = y>(A - A n )> . (1.2) 



n 



By the Banks-Casher relation]] |IJ (taking here the limit V4 — > 00), 

m = ± lim 2m r dA-^L = 2^ , (1.3) 
V±m-*o Jo X 2 + m 2 V4 

this spectral density, evaluated at the origin, is directly related to the appearance of a chiral 
condensate. It follows that the average spacing between eigenvalues becomes roughly constant 
nearA = 0f§, 

AA JL— . (1.4) 

Magnifying p(X) near A = according to the prescription given above thus entails the introduction 
of an associated microscopic spectral density ps(A), which is defined by Q 

PS (X) = lim ±rp(±) . (1.5) 



V4— KX> V4 \ V4 



Applied naively here, ignoring the subtlety of regularization. 
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The precise statement is that this microscopic spectral density should be a universal function. 
The most compelling evidence comes from the fact that the microscopic spectral density defined 
as above, and evaluated in a particular random matrix theory, exactly reproduces the Leutwyler- 
Smilga spectral sum rules [||] and appropriate generalizations ||. 

To understand the significance of this universality conjecture, let us for convenience first focus 
on one example, that of QCD with Nj massless fermions in the fundamental representation. As 
explained in refs. p|, the relevant matrix model partition function is given by a chiral unitary 
ensemble of the form 



Here only symmetry arguments alone have selected the integration to be over complex N x N 
matrices, and the Dirac operator structure is encoded in the determinant. But such symmetry 
arguments alone place no restrictions on the form of ^(A), which appears in the exponent. If 
indeed symmetry arguments alone should determine the microscopic spectral density of QCD from 
this matrix ensemble, then any reasonable choice of V(A) should give the same result for ps(A). 
This universality argument, if correct, allows us to evaluate the integral with a function V(A) 
of our own choice. Some hints supporting this universality had been found in refs. |7|. Clearly, 
the most convenient choice in to make the matrix integral Gaussian except for the determinant 
factor in front. This is the choice made in refs. Q. 

It is important to separate the question of universality of the spectral density of the Dirac opera- 
tor close to the origin into two parts. The first concerns the crucial jump from the d-dimensional 
gauge field integrations to the framework of zero-dimensional matrix models with similar symme- 
tries. The second step concerns the conjecture about the universality of the microscopic spectral 
density within the framework of matrix models. The purpose of this paper is to prove the latter 
of these universality conjectures. The question of whether the universality extends all the way 
from zero-dimensional large- N matrix models to full-fledged quantum field theories such as QCD 
will of course remain unproven. But we shall show, within the context of large- iV matrix models, 
under what conditions the microscopic spectral density and higher correlation functions are uni- 
versal. This is a major step towards understanding the original issue, which concerns the possible 
universality of also the Dirac operator spectrum near A = 0. 

This paper is organized as follows. In section 2 we prove a theorem for the asymptotic behavior 
of generic orthogonal polynomials over a semi-infinite range. By making use of this theorem, we 
compute the associated universal form of the microscopic spectral density as well as all higher 
order correlation functions in the same limit. This universality class is of relevance for four- 
dimensional SU(iV c > 3) gauge theories with fermions in the fundamental representation. In 
section 3 we repeat the analysis for a class of orthogonal polynomials over an infinite range. This 
case is considered to be relevant for three-dimensional SU(-/V C > 3) gauge theories with an even 
number of fermions in the fundamental representation. We conclude in section 4 with a discussion 
about the possible relation between the chiral-flavor (or parity-flavor) symmetry breaking in gauge 




(1.6) 
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theories and the now established universality of random matrices in the microscopic limit. 



2 Chiral unitary ensembles 

In this section we consider the chiral unitary ensemble: 

Z = /dMdet^M exp(-^try(M 2 )l , V(M 2 ) = £ ^M 2k , N f = 0, 1, • • • (2.1) 

where M stands for an (iV + N') x (iV + N') block hermitian matrix whose non-zero components 
are N x N' complex matrices on the off-diagonals, 



M 





W 



(2.2) 



and AM the Haar measure of W. We take N < N' without loss of generality. An ensemble is 
called chiral unitary because of the invariance under the transformation 



W^V^WU, UgV(N), VgV(N'). 



(2.3) 



This model has the same global symmetries as a Euclidean four-dimensional SU(iV c > 3) gauge 
theory coupled to Nf massless fermions in the fundamental representation. The topological charge 
v of the vacuum is identified with \u\ = N' — N (which is kept fixed) and the volume of space-time 
is V4 = N + N' (which is sent to infinity afterwards). Specifically, due to the "U(1)a symmetry" 



{M, 75 } = 0, 75 



In 






Tjv' 



(2.4) 



The partition function ( |2.1[ ) is expressible in terms of the component matrices as well as of the 
eigenvalues after integration over the angular coordinates (U,V) E U(JV) x \](N')/\J(1) N ||, 



all (2N) non-zero eigenvalues M occurs in pairs with opposite signs. 



Z = J dWdet N f (W^W) e 



NtiV(W^W) 



oc 



N 



N 



n( dA * A ?' 



-NV(Xi) 



i=l 



det z 

ij 



det A;." 1 



2(i-l) 



(2.5) 



where we have set a = Nf + \ v\ and suppressed an irrelevant constant of the angular integration. 
The above expression can be interpreted as a positive definite N x N hermitian matrix model in 
H = W^W whose eigenvalues are Ai, • • • , \n > 0. Therefore the problem is reduced to finding a 
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set of orthogonal polynomials Pn a \x) over the semi-infinite interval [0, oo) with the shown weight 
function, 

/■oo 

d AA a e- JVV WpW(A)PW(A) = 5 nm . (2.6) 



o 

The definition of the orthogonality relation (^1]) shows that the problem is well-defined for 
arbitrary real a > — 1. Here we shall only deal with a to be a nonnegative integer, which is 
motivated physically from the above. 

In the orthogonal polynomial method, the polynomials Pn a \\) are usually normalized to be 
monic, 

pH(A) = A™ + -- , (2.7) 

so that the Vandermonde determinant detjj A*" 1 can be substituted by detjj P/"j(Aj) in the inte- 
grand. However, for the present purposes we will employ a different normalization by demanding 
Pl a) (0) = l: 



pW(A) = l + ...+ P i a ¥. (2. 



This is because we shall seek a smooth limiting behavior near A = as n gets as large as N 
(which is taken to infinity), and iV 2 A is kept fixed at the same time (the microscopic limit). We 
thus start by the assumption that Pi°\o) ^ 0, and normalize the polynomials accordingly^ As 
we shall see shortly, this assumption is equivalent to having the origin included in the support of 
the spectral density. 

2.1 Asymptotics of orthogonal polynomials 

We shall prove the following theorem on the asymptotic behavior of these orthogonal polynomials: 
Theorem 1: Let {Pi a) (X)} 

n=o,i,- fre the set of polynomials orthogonal with respect 

to the measure 

dXX a e~ NV ^\ y(A) = ^r^A fc , a = 0,1,..- 

fc>i k 

over the range [0, oo), whose moments are all finite. If the polynomials can be nor- 
malized according to Pb a \o) = 1, then, for fixed x = N 2 X and t = n/N , the following 
limiting relation holds: 



lim Pl a) (^' 



-_ Nt a >(i)vW 1 



where u(t) is determined by 



' We avoid a proliferation of new symbols by employing the S 3,1116 Pji (A) and hn in both conventions 



4 



The proof of this theorem for a = was recently given by one of the authors |]], and is reproduced 
here for completeness. The central idea of how to prove it for arbitrary a = 1,2,... is to use 
induction, starting with the case a = 0. Tracing it back, this means that the Bessel-type behavior 
of the polynomials is a consequence of the three-term recursion relation for the orthogonal poly- 
nomials. As n, N — > oo with t = n/N fixed, this yields a certain "continuum limit" reminiscent 
of the derivation of string equations in the double-scaling limit. The recursion relations thus 
become the defining second order differential equations for Bessel (and Neumann) functions. It is 
the boundary conditions that finally uniquely specify the Bessel solutions. However, for simplicity 
we shall here show the proof by induction. 

We start with the case a = 0. The recursion relation for P n (X) = P?i (A) takes the form (with 
h n = /in 0) and p n = ph )■ 

XP n (X) = -r n (Pn+i(A)-P„(A)-^^i(P n (A)-P n _i(A))} (r n = 

= Y.\n m Pm{\). (2.10) 

m 

This is the well-known three-term recursion relation rewritten in our nonmonic normalization. 
The structure of the coefficients is dictated by a comparison of O(A ) and 0(A n+1 ) terms on both 
sides. The sets of unknowns {h n }, {r n } can be determined iteratively by 



l = -y dA— {e- NV ^P n (X)P n (X)} = NV'CX) nn h n , (2.11) 
= - dA-^ {e- NV ^XP n (X) P n (X)} = (iV(AF'(A)) nn - 2n - l) h n . (2.12) 



In the following we need to know the asymptotic behavior of r n and h n for 

Tl 

n, N — > oo while — = t is kept fixed. (2.13) 



Eqs. Q2.10D , ( |2.11| ) and fl2.12| ) tell us that they should behave as 



/ n\ 1 1 / n \ 1 

r n = r I — j + higher orders in — , h n = — h I —j + higher orders in — . (2.14) 



Then the leading behavior of the matrix A and its powers is 

Xnm — f\ jy) ( ^^nrn $nm+l) ; 

71 v i. x — > / -.of 2/l 



i=-k 



{X k ) nm = r{-f (-)' { k + e ) 6 nm + i (2.15) 
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so that eqs. (|2.11|) and (p.!2|) read 



k-l 



5?» (><'>* 



*(*)' 



Note that eqs. ( [2.16 ) and ( p. 17 ) imply a universal relationship among total derivatives, 



(1/ 2/ (1 ( ^\ + idr = d ( ^ 



(2.16) 
(2.17) 

(2.18) 



Next we expand the right hand side of the recursion relation ( |2.10| ) in terms of 1/iV in the limit 
(|2l^ ) to get 



Ar2l „, A7 x , , . \d 2 P hit) fd r(t)\ dP 
N 2 X P(n, N, A = -r(t) { -—- + -4tt tttt ~t 
y ' ' ; w } dt 2 r(t) \dth(t)J dt 



where the argument N in P(n,N,X) = P n (X) is to indicate explicitly the dependency via the 
coefficient in front of the potential. Only the leading terms up to 0(1/N) have been kept here. 
Eq. ( |2.19| ) tells us that that the arguments of the polynomial appear only in the combinations 
t = n/N and x = N 2 X in the large- N limit. The rescaled eigenvalue coordinate x is to be 
fixed finite hereafter, and is regarded as a parameter in the ordinary differential equation in t. 
Performing the change of variable 



t I— ► u(t) = 



h(t) 



df 



vMO 



(2.20) 



using the relationship ( 2,18; ), and neglecting higher order terms in 1/N, eq. ( 2.1 9| ) reduces to 
Bessel equation of zeroth order: 



1 d d \ „. 

— — u- h x P(u, x) = 0. 

udu du J 



(2.21) 



The unique solution to this equation which satisfies the following boundary condition at t = 
n/N = (u(0) = 0), 

P(0,x) =P (A) = 1, (2.22) 



is 



P(u, x) = Jo (nVx) • 



(2.23) 



We now proceed for a generic integer a by induction. Given a set of polynomials {Pn (A)} 
which are orthogonal with respect to the measure dA A Q e^ NV W on [0,oo), and are normalized 
by Pn a \o) = 1 as in eq. ( |2.g|) , we note that 

p(a+l) (A) _ PitW-Pt\X) ( ^ 24) 
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are polynomials of order n. They are orthogonal to X k , k = 0, • •• ,n — 1 with respect to the 
measure dA A a+1 e~ NV M fu|: 

dXX a+1 e- Nv ^Pi a+1 \X) X k 

dA X a e~ NV{x) (Pi+\(A) - p( a) (X)) X k = 0. (2.25) 

Under the assumption that Pn° +1 ^ (0) 7^ 0, the orthogonal polynomial normalized to P^ +l ^ (0) = 1 
is thus given by 

P^+D(A) = q W . (2.26) 
The microscopic limit of the contiguous relation ( 2.24| ) reads 

p & ^N2 ) = ^¥t N ' { W )+ higher orders in iv" (2 ' 27) 



If we substitute eq. Q2.9| ) for a given a into the above, we confirm that 

5(1+1) / a; a! J a+ i (u(t)yffl 

^(0) = -iV ? (2.29) 
(a + l)V r (*) 

Thus eq. (|2.9| ) holds also for a + 1 after normalizing according to ( |2.26| ) . Together with the result 
( 2.23| ) for a = 0, this completes the proof. 

We have worked out an alternative proof of the same theorem, which does not rely on induction. 
Rather, it is possible to convert the recursion relation for pjf\x) directly into Bessel equation of 
order a, ( ^i+i g^^ 2 " +1 g^ + P^ a \u,x) = 0. Since it is lengthy, we do not reproduce it here. 

2.2 Universal correlations 

Theorem 1 can be used to establish a remarkable universal form of spectral correlators of the 
model (|2.5D. We shall consider the spectral density p(X), and higher order correlators, in the 
same microscopic limit as above. 

To begin, we recall the expression for the integration kernel K^(X,p) associated with the eigen- 
value problem for the positive-definite hermitian matrix H = W^W, 

K N (X,p) = (A/^e-f OWE.)) i. ST P t\X)P^\ti 



i=0 



{Xfi) e ivTw —p '^ Uj 

N—l r 
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Here use has been made of the Christoffel-Darboux identity. In the large-iV limit we may drop 
the superscript a in r n °^ = — Pn Ip^+x- This is because the general recursion relation (the 
counterpart of eq. ( 2.12 ) for general a) 

= N(XV'(X)) nn - 2n - a - 1 (2.31) 

which determines A a \t) leads, in the large- N limit ( |2.13| ), to the same ( p.!7| ) for any finite a. 

On the other hand, the norms of the polynomials can be evaluated iteratively: 



J, 

"n 

We therefore obtain 
1 , 



d\\ a+1 e- NV ^Pt +1) (X)Pt +1) (X) 
We""™ (-PW(A)) ^P^(A) 

■(a+1) 



Ph 



h (a) 
I In 



P 



(a+1) 



(0) 2 



a + 1 V 
Nu(t)/2J 



AT2a+l 



[a 



n 2 



u(t) 



2a 



/i(i) + higher orders in — (h{t) = 



(2.32) 
(2.33) 

(2.34) 



Using Theorem 1 and inserting eq. ( 2.34 ) into ( 2.30 ), we obtain a universal form of the kernel 
(the Bessel kernel of order a) in the microscopic limit 

1.x y u(l) y/xJ a+1 {u(l)y/x)J a {u(l)y/y) - J a (u{l)y/x)y/yJ a+ i(u(l)y/y) 

^N Kn{ N^^ ) = ~2 —y ' 

(2.35) 

The only way this asymptotic kernel depends on the potential V(X) is through the parameter 
u(l). As we shall show now, even this dependence is of a highly universal form. Let us compare 
eqs. ( |2.16D and ( 2.17| ) at t = 1 with the explicit expression for the large-iV spectral density p(z) 
having support on a single interval [—a, a] O]. 



P{z) 
1 



V a 2 — z 2 



2tt 

2k\ fa 2 
k 



k n=0 
k 

= 1 . 



2n\ I a \ 



n 



z 



2fc-2n-2 



This enables us to relate the parameters r(l) and u(l) to a and p(0), respectively, 

a = 2^0) , p(z=0) = 



(2.36) 
(2.37) 

(2.38) 
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The kernel (2.35) therefore only depends on the potential V(A) in the indirect way of setting the 
scale p(0). It has clearly been the assumption throughout that p(A) 7^ 0. This was the starting 
point of the physical motivation from field theory (spontaneously broken chiral symmetries, which 
through the Banks-Casher relation fll.3| ) entails a non- vanishing spectral density at the origin). 
Moreover, it is known that the critical condition p(0) = for a transition where the intervals 
of support move away from the origin is equivalent to Pjv(O) = for the monically normalized 
polynomials If we go back to our proof of Theorem 1, this was precisely the condition we had 
to impose in order to be able to find a smooth continuum limit of the orthogonal polynomials 
in the large-iV limit. Thus, under the assumption that the normalization ( [2.8D is possible, the 



constants u(l) and r(l) are determined to be positive from the identification Q2.38 ). 

It finally remains us to relate the s-point correlation functions of (pn) of eigenvalues of 
H = W^W (or M) to the kernel K N in eq. fl2.35|) . The former is defined as 

a N (\ 1 ,---,\ s ) = (f[±-tr5(\ a -H)) = det K N (\ a ,X b ), (2.39) 
\ a =i N I ^ a ^ s 

p N (zi, ■ ■ ■ , z s ) = / f[ ^tr 5{z a — M)^ = \zi\ ■ ■ ■ \z s \ a N (zf, • • • , z 2 s ), (2.40) 

respectively. Therefore all the formulae for their microscopic limits (x a = N 2 X a , ( a = Nz a fixed) 
**(*!, • • • , x.) = lu^ jj;'N(%5, ■ ■ ■ , (2-41) 
Ps(Cu ■■■Xs)= to^PNijj, •••'§) = ICil • • • IC.I ^s(Ci 2 , • • • , C 2 ), (2.42) 

that were previously calculated for the Laguerre (in the -ff-picture) or chiral Gaussian (in the 
M-picture) unitary ensemble [^, |i"2"| , hold universally. 

In particular, the spectral density of the chiral unitary ensemble 

Pn(z) = (±tr6(z-M)} = \z\K N (z 2 ,z 2 ) (2.43) 

takes the universal form 

Ps(C) = MO)) 2 ICI (^(2^(0)0 - J a+ i(27rp(0)C)Ja-i(27rp(0)C)) (2.44) 

in the microscopic limit. As expected from the Gaussian case, a matching condition between the 
microscopic and macroscopic (ordinary large- iV) spectral densities is satisfied: 

lim p s (()=p(0). (2.45) 
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3 Unitary ensembles 

In this section we consider the unitary ensemble: 

Z= [dMdet 2a Me- m ^ M2 K V(M 2 ) = ^ ^-M 2k , a = 0,l,--- (3.1) 



fc>i 



where M stands for an N x N hermitian matrix and dM its Haar measure. This model shares the 
same global symmetries with a Euclidean three-dimensional SU(iV c >3) gauge theory coupled to 
Nf = 2a massless fermions in the fundamental representation The partition function fl3.1|) 
is expressible in terms of the eigenvalues after integration over the angular U(iV) coordinates 



Z= fnfdA, X 2a e~ NV ^ 

J -co ,. =1 



detAr 1 



'•j 



2 



(3.2) 



Therefore the problem is reduced to finding a set of orthogonal polynomials Pn (A) over the 
infinite interval (— oo, oo) with the shown weight function, 

dAA 2a e -W(A2)pM (A) p(a) (A) = h («) ^ (3.3) 

3.1 Asymptotics of orthogonal polynomials 

We shall prove the following theorem^ on the asymptotic behavior of these orthogonal polynomi- 
als: 

Theorem 2: Let {Pi a) (X)} 

n=o.l,- be the set of polynomials orthogonal with respect 

to the measure 

dXX 2a e- NV ^ 2 \ y(A 2 ) = £^A 2fc , a = 0,l,--- 

fc>i 2k 

over the range (—00,00), whose moments are all finite. If the polynomials can be 
normalized according to P^„ (0) = P^n+i'ity = 1; then, for fixed x = NX and t = 
2n/N , the following limiting relations hold: 



lim 



X 



lim NP£U—) 



1 J a -l (u(t)x) 
= T(a + -) 2 r , (3.4a) 

n=Nt/2 2 (u(t) x /2) a -2 

3 J , 1 (u(t)x) 
= x T(a + -) + 2 (3.4b) 

n=Nt/2 2 ( u (t) x /2) a+ 2 



where u(t) is determined by 



* Theorem 2 can be regarded as an extension of Theorem 1 for half-integer a. 
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The proof of this theorem for a = was sketched in ref. fL4|| , and below we shall elaborate it in 
a rigorous form. The recursion relation for monically normalized polynomial 

P n (A)^p(°)(A) = A n + --- (3.5) 

takes the form 

AP n (A) = P n+1 (A) + 

fn -Pi— 1 

(A), (3.6) 

with 

r n = . (3.7) 



Since the monic normalization ( |3.5|) does not distinguish the parity of n, neither does r n . Eq. fl3.6|) 
immediately implies 

= P 2 „+ 2 (0) + r 2n +i P 2n (0). (3.8) 

Due to the Z 2 symmetry of the weight function, the polynomials have definite parities. Therefore 
we are lead to change the normalization of polynomials as follows: 

Then the recursion relation ( |3.6| ) in this normalization, applied twice to P 2n (A), takes the form 

A 2 P 2 „(A) = -r 2 „+i(P 2n+2 (A)-P 2n (A))+r 2 „(P 2n (A)- J P 2n _ 2 (A)) 

= E(^)2",2 m P 2m (A). (3.10) 

m 

Here use is made of eq. (|3.8j). A similar equation for P 2n +i(A) expresses the matrix elements 
(A 2 ) 2n +i,2m+i hi terms of the r n 's. The coefficients {r n } are iteratively determined by 

° = -7^ I" dX^-{e- NV ^XP n (X)P n (X)\ = 2N(X 2 V'CX 2 )) nn -2n-l. (3.11) 
n n J-oo dX L J 

The asymptotic behavior of r n is determined in the large- limit, 

2n 

n, N — » oo while — = t is kept fixed, (3-12) 

r n = r (yj^j + higher orders in — . (3.13) 

\T,9 k ^)r{tf = t. (3.14) 



to be 



Here r(f) is given by 



Next we expand the right hand side of the recursion relation ( [3. 10 ) in terms of 1/N to get 
(P(2n,N,X) = P 2n (X)) 

N 2 X 2 P(2n,N,X) = - (m*)^ + 2^|^) = ~ (2 v /^)^) 2 J P(2n,iV,A). (3.15) 
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It tells us that that the arguments of P appear only in the combinations t = 2n/N and x = NX 
in the large- N limit. After the change of variables 



1 1 ^ u(t) = [ 
Jo 



At' 



we equivalently have 



du 2 



o 2^H¥)' 



+ x 2 P(u,x) = 0. 



(3.16) 



(3.17) 



The unique solution to this trigonometric differential equation which satisfies the boundary and 
parity conditions 

P(0,x) = l, P(u,-x) = P{u,x) (3.18) 



is 



P(u,x)= lim P 2n (-) 



n=Nt/2 



1 J-i (ux) 

COSUX = T(-) j- 

2 (uxjiy^ 



(3.19) 



The odd-order polynomials are constructed out of normalized even-order ones as follows. 

-P2n+2(A) — -P2n(A) 



P2n+lW 



A 



(3.20) 



are odd polynomials of order 2n + 1, and are orthogonal to X 2k+1 , k = 0, • • ■ , n — 1: 

dAe-^Wp 2n+1 (A)A 2fc+1 



-oo 
oo 



dAe -w(A) (P 2n+2 ( A )-P 2n (A))A 



2k 



0. 



The microscopic limit of ( |3.20| ) reads 

~ x 
^2n+l(jy) 



- 29 P ( X ) 

n=Nt/2 Xdt N 



n=Nt/2 



(3.21) 



(3.22) 



Under the assumption that P2n+i'(0) / 0, change of normalization to P2n+i'(0) = 1 is imple- 
mented by 

-FWt+i (A) 



2ra+l 



(A) 



(3.23) 



iWi'(o) 

By inserting the asymptotic form of even-order polynomials ( p. 19 ) into (3.22) and normalizing 
according to ( |3.23 ), we conclude 



lim NP 2n+1 (^- 
N— >oo iv 



smux „,3 N J\ {ux) 

= = £r(-j — r 

n=Nt/2 U 2 ( ux /2)2 



(3.24) 
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Now we proceed for a generic a by induction. Given a set of polynomials {P&*\\)} which are 
orthogonal with respect to the measure d\ \ 2a e~ NV ^ x ) on [—00,00], and are normalized by 

P^(0) = P 2 ( n+i'(°) = we note that 



A 2 



(3.25) 



are polynomials of order n. They are orthogonal to each other with respect to the measure 
dAA 2(a+i) e -W(A 2 )_ Unc ier the assumption that P 2 ( n +1) (°)> P2n+l'{ ) + °> the correctly normal- 
ized orthogonal polynomials are given by 



it +1) (A) 



P 



(«+i) 



2d 



(A) 



P 



(a+l) 



P 



( a + 1 )(c^\ , 2n+1 



(A) 



P 



(o+l) 



2n+l 



(A) 



2n '(0)' ~" ' ' P£J?'(0) 

If we substitute eqs. (3.4) for a given a into the microscopic limit of eq. ( |3.25| ), we confirm that 

T(a + \ 



(3.26) 



p(«+i) r x > 



p(a+l)/^ 



-JV- 



Ja+| («(*)«) 



3(«+l) 

2n 



P 



(o+l)/ 
2n+l 



(0) 



(0) 



n=Nt/2 

n=Nt/2 
n=Nt/2 
n=Nt/2 



v 7 ^) (u(t)/2) a -? x a+ ^ 



T(q + |) Jg+I i.u(t)x) 
VW) (u(t)x/2) 



a+l 



-N 
-N 



u(t) 



(2a + l)v / Kt) 
u(t) 



(3.27a) 

(3.27b) 
(3.28a) 
(3.28b) 



(2a + 3)v^*J' 

Thus eqs. (3.4) hold also for a + l after normalizing according to (|3.26j ). Together with the results 
( |3.19| ) and ( |3.24| ) for a = 0, this completes the proof of Theorem 2. 



3.2 Universal correlations 

The integration kernel associated with the eigenvalue problem of M is (for brevity we set ./V 
odd = 2Af + 1) 

-, N—l d (q) 1 \ \ D (a) / \ 



/,(«) 

n=0 'm 
AT \ D (a) (\\-o( a ){,,\ r>( a ) ( \\ r>( a ) 



= (AM) a e-*^HV0H))^j5- 



1 \^i(a)^(m)-^a7(A)p 2 v a7 +1 (m) 



vn=0 " 2 n 



(a) 



A — /i 



(.3.29) 



Here use is made of the Christoffel-Darboux identity, and the proportionality constant in the 
second line is determined by matching O(A ^°) terms. Noting that the orthogonality relation 
in the microscopic limit 



N 2c 



n=Nt/2 



-,(«), x_- 
2m \ jy ' 



m=Nt'/2 



4? pit - 



(3.30) 
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is identical to the Bessel closure relation (inversion of Hankel transform) 



1 

dxxJ i(ux)J i(ux) = —5(u — u), (3.31) 

2 2 U 



we can determine the asymptotic form of the norm as 



h 2n = J^ T ( a + 2> V r ^ + hi S her orders in AT (3 ' 32) 



Accordingly we obtain 

^ J_ _ iV 2a JV f 1 u{t) 2a N 2a+1 u(l) 2a+1 

_ r(a + i) 2 2 2 -' 2 io t 7?t!j"r(a + i)2 22« 2a + 1 ' (3 ' 33) 



Using Theorem 2 and inserting eq. ( |3.33j ) into ( |3.29|) , we obtain the universal form of the kernel 
(the generalized sine kernel of order a) in the microscopic limit 

1 x y u (l) J a+ i(u(l)x)J a i(u(l)y)- J a i(u(l)x)J a+ i(u(l)y) 

jfe, n Kn ^n ) = J T^^ 2 x~^ • (3 ' 34) 

Therefore all formulae for the correlation functions 

p N (X u ---,X s ) = ( f[ i-tr<5(A a -M)\ = det K N (X a ,X b ), (3.35) 

\ —i I l<a,o<s 



previously calculated for the Gaussian unitary ensemble in the microscopic limit [15, 12] 

Ps (xi, • • • , x s ) = lirn^ WsP n(jj, ■ ■ ■ , jj), (3.36) 

hold universally. Specifically, the spectral density Pn(X) = {1/N tr8(X — M)) = Kjy(X,X) univer- 
sally takes the form (now p(0) = u(l)/ir) 

Ps(x) = (^) 2 x (j 2 a+l + Jl_ h - J a+ iJ a „l ~ J a -\J a+ l) WOW • (3-37) 
It enjoys the matching condition 

lim ps(x) =p(0). (3.38) 



4 Conclusion and speculation 

In the present work we have proven two theorems concerning the asymptotic behavior of orthog- 
onal polynomials for a broad class of measures. They show a remarkable asymptotic universality 
near the origin, a relation previously only known for the classical orthogonal polynomials of gen- 
eralized Laguerre and Jacobi type |l(J. Furthermore, we have used this highly universal behavior 



14 



to prove universality conjectures regarding the spectral densities and their s-point correlators in 
the microscopic limit. 

All of these proofs refer to the zero-dimensional matrix model language alone. And, as we have 
emphasized in the Introduction, the assumption that this universality can be applied to full d- 
dimensional quantum gauge theories with fermionic fields is a priori very far from obvious. Of 
course, our matrix model universality theorems are interesting in their own right. But from a 
particle physics perspective the real focus should be on the application of these theorems to the 
Dirac operator spectrum Q. 

How could we imagine the link between random matrix models and gauge theories established 
in this connection? Consider a Euclidean four-dimensional SU(A C ) gauge theory coupled to Nf 
light fermions in, say, the fundamental representation. Assume that Nf is small enough to allow 
for the conventional scenario of spontaneous chiral symmetry breaking. Let us simply call such 
a theory "QCD" . Its partition function can be written formally as 

Nf 

Z QC d = I 77 f^ T LI WfWf e-W-J*t&»"ty*f , (4.1) 
J (Gauge) ^ 

where S[A] is the gauge field part of the action. Consider the theory in a finite volume V4. Let 
us now seek a low-energy effective description of this theory. One way is the chiral Lagrangian 
(xL) approach, which in this case introduces considerable simplification due to the suppression 
of all derivative terms (since one is interested only in the very soft modes) ||]. Here the chiral 
condensate at zero mass, {ipip} = £ is by construction assumed to be non-vanishing. Below is a 
schematic picture of how one could imagine random matrix theory to fit into this framework: 

Z QCD = I (cg^y e~ SlA] Uf det(ip+ m f ) H Z xUE = / dMe~ NV ^ Uf det(M + m f ) 
I (a) (c) I 

Z X L = I sv(Nf) dUe N ^ mUelWJ +^ H Z xGlJE = fdMe- N ^* M2 Ufdet(M + m f ) 



This link (a) is, however, based on the global symmetry breaking pattern SU(Nf)i, x SU(iVf )r — ► 
SU(iVy)v alone and thus by no means refers to the microscopic theory. 

Let us now consider a different route (b)-(c)-(d): suppose we start with QCD in a finite volume, 
and with a sharp ultraviolet cut-off. (One would like to think of conventional lattice regulariza- 
tions, were it not for the known difficulties of defining massless fermionic degrees of freedom in 
that case.) Consider the formal change of integration variables from the gauge field A fl (x) to 
the Dirac operator \If) = A(x) + if) = M itself (b). Since V4 is the size of the matrix M, we 
write N = V4. Although we can hardly imagine computing the Jacobian J of this transformation 
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the resulting partition function would formally look like the random matrix model above, with 
V(M) = S[A(M)] — log J{M). Once phrased in a matrix model language, we can apply our Theo- 
rem 1 to substitute the uncomputable V(M) by another simple measure, say that of the Gaussian 
unitary ensemble (GUE) with V(M) = S 2 trM 2 . This GUE can be shown to be expressible as 
the zero-dimensional chiral Lagrangian in the microscopic limit (d) Q]. If the relation (b) be- 
tween a properly regularized version of QCD and a particular random matrix model were not 
still very formal and weakly established, the diagram would otherwise close itself full circle. The 
diagram above at least serves as an intuitive picture outlining the logical routes between the dif- 
ferent formulations, and the role played by our universality theorems proposed in this connection. 
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